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ABSTRACT 

We report on the first 180 days of RXTE observations of the outburst of the black hole candidate 
ICR J17091-3624. This source exhibits a broad variety of complex light curve patterns including 
periods of strong flares alternating with quiet intervals. Similar patterns in the X-ray light curves 
have been seen in the (up t o now) unique black hole system GRS 1915-1-105. In the context of the 
variability classes defined bv iBelloni etall (|2000f ) for GRS 1915+105, wc find that ICR J17091-3624 
shows the v, p, a, A, l3 and /i classes as well as quiet periods which resemble the x class, all occurring 
at 2-60 keV count rate levels which can be 10-50 times lower than observed in GRS 1915-1-105. The 
so-called p class "heartbeats" occur as fast as every few seconds and as slow as ~1GG seconds, tracing 
a loop in the hardness-intensity diagram which resembles that previously seen in GRS 1915-1-105. 
However, while GRS 1915-f 105 traverses this loop clockwise, ICR J17091-3624 does so in the opposite 
sense. We briefly discuss our findings in the context of the models proposed for GRS 1915+105 and 
find that either all models requiring near Eddington luminosities for GRS 1915+105-like variability 
fail, or IGR J17091-3624 lies at a distance well in excess of 20 kpc or, it harbors one of the least 
massive black holes known (< 3Mq). 

Subject headings: X-rays: binaries — binaries: close — stars: individual (IGR J17091-3624, GRS 
1915+105) — Black hole physics 



1. INTRODUCTION 

Observations with the Rossi X-ray Timing Explorer 
(RXTE) have led to extraordinary progress in the knowl- 
edge of the variability properties of many different types 
of sources, particularly of black hole candidates (BHCs) 
and neutron stars (NSs) in low-mass X-ray binaries (e.g, 
Ivan der Klij |2006[) . Both BHCs and NS are known 
to exhibit distinct "accretion states" (usually defined 
in terms of their X-ray spectral shape and variability), 
whose characteristics are thought to be intimately re- 
lated to the physics of the accretion flow and its in- 
terac t ion with the c entral compact o bject (e.g.. iBellonil 
20101: [Rmfiilard fc McChntock 200& Ivan der Klisli2Q0a 
Belloni et al.ll2011l) . 

Of all Galactic BHC X-ray binaries known today, 
one of the most prolific in terms of state transitions is 
GRS 1915+10 5. In outburst since its discovery in 1992 
with WATCH (jCastro-Tirado et al.|[T992,). it is a 33 days 
orbital period binarv svstem jGreiner et al.l l200l!) har- 
boring a 14 ± 4.4 black hole (jCreiner et al.ll200ll: 

tirlaftis & Grcincr 2 0(M). At a distance of ~ 12.5 kpc 
lirabel fc R,odrfgue3TT99l . GRS 1915+105 is very of- 
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ten at Eddingto n or super-Eddington luminosity (e.g., 
iDone et al.ll2004D . 

GRS 1915+105 is known to show quasi-periodic oscil- 
lations (QPOs) in the 0.001-70 Hz range, some of which 
are t he same as those generally seen in other BHCs 
(Mor gan et "all 119971: iMarkwardt et"an]1999t iReig et al.l 
iSoOOi : IStrohmaveill2001l : ISoleri et al.ll2008D . However, so 
far GRS 1915+105 has been unique in that its X-ray 
light curves exhibit more than a dozen different pat- 
terns of variability usually called "classes" (which are re- 
ferred to with Greek letters) , most o f which are high ani - 
plitude and highly-structured fe.g.. IBelloni et alll2000D . 
Most of this structured variability is thought to be due 
to limit cycles of accretion a nd ejection in an unstable 
disk (see, e.g IBelloni et al.llT997: Mirabel et al.l 119981: 
iTagger et all 12004 iNeilsen et al.l 1201 iL and references 
within). Much effort has been expended in order to un- 
derst and why GRS 1915+10 5 is so unusual among BHCs 
(e.g. jFender fc Belloni|[200^ . It has been proposed that 
the high accretion rate estimated f or GRS 191 5+105 
might be the determining factor (e.g.. lDone et al.i r2004): 
however, the lack of any other source showing similar 
characteristics has prevented definite conclusions. 

ICR J 1709 1-3624 was discovered with INTE- 
GRAL/IBIS durin g a Gal actic Center observation on 
April 2003 (Kuulkc rset al.l l2003:i and later again in 2007 
()Capitanio et al.. .20091 ). Reexamination of archival data 
from different missions showed that I CR J17091-3624 
was al so active in 199 4 , 199 6, 2001 (Rcvnivtscv et alj 
120031 lin't Zand et"al] l2003t ICapitanio ct al. 2006). 
Although a neutron star could not be excluded based 
on spectral ch aracteristics and the o utburst radio/X- 
ray flux ratio, ICapitanio et al.l ()2006f ) concluded that 
IGR J17091-3624 is most probably a black hole. A new 
outburst was detecte d with Swift/BAT in February 2011 
(|Krimm et al.l 1201 ll ): the radio/X-ray characteristics 
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Fig. 1. — Upper Panel: Swift/BAT light curve. Lower Panel: Circles: RXTE Crab normalized 2-20 keV intensity. RXTE values are 
averaged per observation, background subtracted and dead-time corrected. Triangles: RXTE 2-50 keV flux RXTE observations between 
day 55601 and 55615 are affected by the nearby source GX 349-1-2 (Sco X-2); the intensity during this period was arbitrarily fixed to 65 
mCrab and no X-ray flux is reported. Shaded areas mark quiet periods (see text). 

during the first ^-^40 days, combine d with the discovery 
of QPOs (jRodriguez et al.l |20TI3), further suggested 
that IGR J17091-3624 is a black hole. 

Recently we reported the discovery of 10 mHz 
QPOs in RXTE observ ations of IGR J17091-3624 
(|Altamirano et ahl l2011aD very similar to those in 
GRS 1915-1-105. The suggested link between these 



sources was strengthened by our discovery of (i) a contin- 
uous progression of regular, quasi-periodic flares occur- 
ring at a rate of 25-30 mHz (jAltamirano et al.ll2011bD . 
and (ii) a broad vari ety of complex patterns a lternating 
with quiet intervals (jAltamirano et al.]|2011bl B). resem- 
bling respectively, the so called p class ( "heartbeat" ) os- 
cillations and the complex /3 class patterns, both so far 
seen only in GRS 1915+105. 

The existence of a source showing similar X-ray vari- 
ability to that seen in GRS 1915-1-105 opens a new win- 
dow of opportunities to understand the physical mecha- 
nism that produces the highly structured X-ray variabil- 
ity. Multi-wavelength and high spectral resolution stud- 
ies of IGR J17091-3624 as compared with GRS 1915-^105 
can help gaining further insights into the role o f disc-jet 
coupling (e.g., review by 'Fender &: Bellonil l2004) and ac- 
cretion disk winds (e.g, Ncilscn ct al. 201l|) in accreting 
BHCs. Furthermore, it allows the possibility to test the 
role of the accretion disk size, evolutionary state of the 
companion star, and the evolution of the disc structure 
as an explanation of the long- and short-term X-ray vari- 



ability seen in GRS 1915+105 fe.g.. lDone eraIll2004D . 

In this Letter we use RXTE data to describe the 
phenonienological similarities between IGR J17091-3624 
and GRS 1915+105. We demonstrate that IGR J17091- 
3624 shows the same type of high amplitude and 
highly-structured variability previously observed only in 
GRS 1915+105, although at much lower count rates. The 
analysis presented here is the first step of a larger pro- 
gram to compare IGR J17091-3624 and GRS 1915+105 
in detail using data obtained with RXTE, Swift, XMM- 
Newton as well as optical facilities. 

2. OBSERVATIONS, DATA ANALYSIS 

IGR J 1709 1-3624 is observed with the Proportional 
Coun ter Array (PGA; IZhang et aLlll993t iJahoda et al] 
|200(J ) on-board RXTE almost daily since the outburst 
began in February 2011 (at the time of submission of this 
Letter, IGR J17091-3624 was still active). We use the 
first 147 observations, covering ^180 days. We also used 
1700 archival RXTE observations of GRS 1915+105. 

Power spectra and light curves were p roduced from 
the PGA using standard tech niques (e.g., iBelloni et al.l 
120001 : lAltamirano et al.|[2008h . Deadtime corrected en- 
ergy spectra averaged over each observation were cre- 
ated from PCU 2 data; a systematic error of 1% was 
added to all channels. Response matrices were created 
with PCARSP (V11.7 . 1) us ing the position reported by 
iKennea fc Capitaniol ()2007f) , thereby taking into account 
a constant 0.4° pointing offset used after February 23rd, 
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Fig. 2. — The upper and lower panels of each frame show a light 
curve of IGR J17091-3624 and GRS 1915+105, respectively. Count 
rates are given in 1 s bins, per PCU, 2-60 keV and background sub- 
tracted. The IGR J17091-3624 A-C intervals last Tl seconds and 
come from observations 96420-01-05-00, -06-00 and -07-01, respec- 
tively. For GRS 1915-f 105 they last T2 seconds and come from ob- 
servations 10408-01-40-00, 20402-01-34-00 and 93701-01-02-01, re- 
spectively. R = MINgrs 1915+105 /M I Njgr ji7091~3624, where 
MIN is an average count rate during minima in the above light 
curve and is given for each frame as an order of magnitude estimate 
of the flux ratio. Greek letters on the side indicate variability the 
class. 
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Fig. 3.— Similar to Figure [2] The IGR J17091-3624 D-F in- 
tervals come from observations 96420-04-03, -08-03 and -09-06, re- 
spectively. For GRS 1915-1-105 they come from observations 20187- 
02-01-00, 95701-01-31-00 and 10258-01-10-00, respectively. 
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Fig. 4. — Left panels show the hardness-intensity diagram for flares observed during the p variabihty class in IGR J17091-3624 (top; 
ObsID 96420-01-04-03) and GRS 1915-1-105 (bottom; ObsID 96378-01-01-00) occurring at an average period of Tl=70.96 seconds and 
T2=63.72 seconds, respectively. Arrows mark the time evolution. Inset shows representative flares. Light curves and colors are estimated 
from 1 sec averages. Intensity is the count rate in the 2-60 keV range (absolute channels 0-240) and hard color is the 6.5-15.0 keV / 
2-6.5 keV count rate ratio (channels 15-35 and 0-14, respectively). Right panels show representative power spectra from averages of 512 sec 
segments during the p variability class for IGR J17091-364 (top, ObsID:96420-01-05-000, MJD 55647.9) and for GRS 1915-1-105 (bottom, 
ObsID:40703-01-07-00, MJD 51235.3). 

2011 (MJD 55615) to avoid a bright nearby source (see 
below). All spectra were satisfactorily fitted in the 3-25 
keV band with an absorbed disk-blackbody plus power 
law model ( phabs* (diskbb + powerlaw)) using Xspec 
(lArnaudl[l99 6) 



12 (lArnaudI 1199 6) with N,, fixed to 1.1 x W cm"^ 
(|Krimni et al.ll201l[ ): some spectra required a Gaussian 
at 6.5 keV for a good fit. 

3. RESULTS 

Figure [T] (upper p anel) shows the Swift/BAT 
(|Barthelmv et al.l |2005() daily light curve of the first 
200 days of the outburst (assuming outburst onset on 
MJD ^^55580). The outburst reached a maximum of 
0.026 counts/cm^/s (15-50 keV) ~27 days after onset; 
then decreased and remained below 0.01 counts/cm^/s. 
The lower panel of Figure [1] shows the 2-20 keV, Crab- 
normalized intensity and the 2-50 keV un-absorbed flux 
from each RXTE/PCA pointed observations. RXTE 
observations started 21 days after onset. The first 
10 observations (until MJD 55615) are affected by the 
contribution of the bright, variable source GX 349-1-2 



(Sco X-2) which was w ithin the 1° PGA FoV (see, also, 
[Rodriguez et alj [201131) . For reference, in Figure [1] we 
show these observations at an arbitrarily fixed intensity 
of 65 mGrab; we do not report their X-ray flux. 

From the variability point of view, we find that dur- 
ing MJD 55601-55632 , the X-ray light curves were well 
characterized by broad band noise and a QPO moving 
in time from 0.1 Hz up to '^6 Hz. The evolution of 
the power spe ctral components (see, iPahari et al.l 120111 : 
IShaposhniko^l2011l : [Rodriguez et al.ll2011al for more de- 
tails) is similar t o that seen in BHG transients in the 
hard state (e.g., lBelloni|[2010D . further suggesting the 
black hole nature of IGR J17091-3624. 

On MJD 55634 IGR J17091-3624 first showed a ~10 
mHz QPO. Since then, a broad variety of complex pat- 
terns occurred, including strong variab ility in the form of 
flares alternating with quiet intervals (jAltamirano et aP 
l2011bl lH). The variability shows re markable similarities 
with that of GRS 1915-M05 (e.g., iBelloni et aLll2000D 
suggesting a common mechanism. In Figures [2] and [3] 
we show six different examples; top and bottom frames 
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of each panel correspond to data from IGR J17091-3624 
and GRS 1915+105, respectively. 

Based on the classification of iBelloni et al.l ()2000D . 
panel A shows segments of the v class, and panels B 
and C two different varieties of p class. The v and p 
variability classes in GRS 1915+105 are characterized by 
quasi-peri odic "flares" recurring on time scales bet ween 
~40 fe.g.. iMassaro et al.l 120101: iNeilsen et al.ll20lH) and 
~120 seconds (jBelloni et al.ll2000l ). The main differences 
between the v and p classes are that (i) the shape and 
period of the flares in v can be more irregular than in p 
and (ii) that the v flares show a characteristic structure 
in their profiles, notably a dip followed by a secondary 
peak after the main one (see Figure [5]). 

Panel D shows segments of the a class, in which 
"rounded-bumps" are sometimes accompanied by sharp 
peaks which last a few seconds. Note that in our ob- 
servations of IGR J17091-3624 we do not see the typical 
'^1000 s quiet interval which precedes the bumps in the a 
class of GRS 1915+10 5 (not shown, see, e.g.. Figure 19d 
in Belloni et al.ir2000l) . Panel E shows light curves from 
either the /3 or the A class. These light curves consist 
of quasi-periodic alternation of low-quiet periods with 
highly- variable and oscillating ones. Panel F shows seg- 
ments of the p class, where we observe periods of very 
rapid (factor of 2-4 changes in count rate in less than 5 
seconds) and almost incoherent variability. 

In addition to the light curves shown in Figures [5] and [3] 
, IGR J17091-3624 displays complex light curves com- 
bining characteristics of different classes. It also shows 
quiet periods in which 1-scc light curves are flat (< 10% 
fractional rms amplitude at < 1 Hz). These periods are 
marked with grey shadowed areas in Figure [T] and coin- 
cide with the times when IGR J 1709 1-3624 is detected 
most strongly by Swift/BAT, i.e. when the spectrum is 
hard. In terms of the GRS 1915+105 variability classes, 
these periods could correspond to the x class. 

Figures [5] and [3] not only shows the remarkable simi- 
larities in the shapes of light curves from both sources, 
it also shows two clear differences: (i) the time scales 
can be different (IGR J17091~3624 tends to be faster, 
see below) and (ii) the average count rate (or flux) 
of the source can be much higher (factor 10-50) in 
GRS 1915+105. Some of the spectral and timing anal- 
ysis performed for GR S 1915 +105 on ti me scales of sec- 



onds (e.g., iMar kward t et al.l 
l.ll20f ' 



iDone etaU mm . ISoleri et al.l 



1999; Bell oni et al.l \200dt 
20081 ) is prevented by the 



combination of relatively low count rate and faster vari- 
ability in IGR J17091-3624. 

Given that we find the v and p variability classes (pan- 
els A-C) in about 35% of our observations (another ^-^35% 
of the light curves are flat, and the remaining ^30% are 
a mix of the a, {3, p, X and unclassified ones) and that 
p is one of the best studied classes in GRS 1915+105, in 
the rest of this Letter we constrain ourselves on further 
comparison of the ly and p classes between sources. More 
detailed comparison of the other variability classes will 
be presented in upcoming papers. 

Some of the observations clearly show only the p or 
the p types of flares, some show a mix and sometimes 
differentiating the two classes is difficult due to the low 
statistics. In any case, these flares can occur as fast 
as every few (2-5) seconds (e.g., ObsID: 96420-01-22-04, 
MJD 55768), and as slow as every ~100 seconds (e.g., 



ObsID 96420-01-03-01, MJD 55634). This means that 
the oscillations in IGR J17091-3624 can be faster than 
those in GRS 1915+105, but not as slow. The fractional 
rms amplitude of the flares covers a range from ^2% to 
up to 50-60%. If one assumes that the minimum pe- 
riod that a quasi-periodic feature can reach scales pro- 
portional to s ome power of the mass of the compac t ob- 
ject (see, e.g.. lBelloni et al.lll997tlFrank et al.ll2002[ ). and 
that a 2-5 s recurrence time of the p/v flares IGR J17091- 
3624 versus --40 s in GRS 1915+105 is due to a differ- 
ence in mass, then our results suggest the black hole in 
IGR J17091-3624 could be a factor of a few less massive 
than the 14 ± 4.4 Mg of GRS 1915+105. 

The variability classes are known to exh i bit d istinc- 
tive spectral evolution (e.g., IBelloni et al.l |2P0tf). In 
a color-color diagram (CD) or hardness-intensity dia- 
gram (HID), one observes that each flare from the p 
(and sometimes i^) class traces a loop (or "ring" , see 
iVilki fc NevalainenI [19981 : iBefloni et aLll2000[ ). Figure H 
(bottom-left) shows the HID for ten consecutive flares 
from a single observation of GRS 1915+105 (inset shows 
a representative flare). The loop in the HID is always 
traversed clockwise. 

The HIDs and CDs from single flares of IGR J17091- 
3624 are dominated by low statistics, so we could nei- 
ther exclude nor confirm loop-like patterns similar to 
those observed in GRS 1915+105. To improve statistics, 
we used intervals where more than 10 fiares occurred 
approximately periodically and folded each interval at 
the best average period. Although this process washed 
out some of the structure in the light curves (due to 
the quasi-periodicity of the signal and the profile differ- 
ences between flares), all the HIDs we created using the 
same energy bands as in GRS 1915+105 (which are fixed 
by the observing modes used) always showed a loop re- 
sembling that seen in GRS 1915+105. In the upper-left 
panel of Figure |4| we show a representative example (in- 
set shows the average profile of the fiare). However, in 
all cases the loop is traversed in an anti-clockwise sense, 
i.e., opposite to what we see in GRS 1915+105. As the 
hardness ratios are a crude characterization of the spec- 
trum, detailed spe ctral modeling of these loops (e.g., 
INeilsen et all I2OIIII are needed to understand whether 
there is a physical difference between the p class seen in 
both sources. 

Figure |4] (right panels) shows representative power 
spectra of the p variability class in IGR J17091-3624 and 
GRS 1915+105 . The power spectra share the same main 
features: (i) a low-frequency QPO (with high harmonic 
content) due to the "flares" and (ii) a QP O with char- 
acter istic frequency between 6-10 Hz (e.g., iMuno et al.l 
Il999f ). In addition, for both sources we sometimes find 
a "bump" with characteristic frequency between 1 and 5 
Hz. 
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4. DISCUSSION 

literature exists attempting 
complex variability observed in 
Some authors propose that the 
(close to, or super-Eddington 
1915+105 is the determining factor 
I1997I: IVilhu fc Ncvalaineg 
Navakshin et al. 2000; .Janiuk et all 
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Fig. 5. — Mass of the compact object versus distance to the bi- 
nary system assuming IGR J17091-3624 is emitting at Eddington. 
The three curves correspond to different (1, 2 and 3) correction 
factors between the 2-50 keV and the bolometric flux. 

ences therein). Although it is generahy accepted that 
the complex X-ray variability in GRS 1915-1-105 results 
from disk instabilities, the exact nature of the instability 
remains unknown. The lack of at least a second source 
showing similar characteristics has prevented definite 
conclusions. 

In this Letter we show for the first time that another 
source, the BHC IGR J17091-3624, can show the same 
broad variety of complex light curves as GRS 1915+105 
(at least 7 of the 12 variability classes observed in 
GRS 1915-1-105, in addition to unclassified ones; note 
that at the time of submission of this Letter IGR J17091- 
3624 is still active) . Although the comparisons presented 
in this paper constitute only a first step, the observed 
similarities suggest that the complex light curves of the 
two sources are produced by the same physical mecha- 
nisms. If true, the low flux of IGR J17091-3624 compared 
with GRS 1915-1-105 combined with the circumstance 
that currently neither the distance to IGR J17091-3624 
nor the mass of its compact object are known, raises 
the fundamental question: is IGR J17091-3624 close to 
Eddington or not at times when showing the same char- 
acteristic X-ray variability? 

A scenario in which IGR J 1709 1-3624 is not emit- 
ting at close to the Eddington but only at a few per- 
cent, is at variance with models where the variability 
is explained as due to disk instabilities that on ly oc- 
cur at high luminosity (e. g.. IVil hu fc Ncy alaincn j 1 19981: 
Belloiii et a l. 2000; Navaks hin et all feoOOt iJaniuk et al.l 
20021: iDone et al.. ,2004; ,Neilsen et all 1201 iL and refer- 
ences therein). However, in this scenario IGR J17091- 
3624 would follow radio/X-ray correlation of some or 
most BHCs (depending if the distance is closer to ~ 11 
kpc or ~ 17 kpc, respectivelv: [Rodriguez et al1l2011b| l. 



A scenario in which IGR J 1709 1-3624 is emitting at 
close to the Eddington limit puts constraints on its mass 
and distance. Although this scenario would imply that, 
independently of the distance, IGR J17091-3624 does 
not follow the sta ndard Radio/X -ray luminosity plane 
for BHCs (Rodriguez et al.ll2011bf ). this would be similar 
to GRS 1915-1-105, as it does not follow the radio/X-ray 
flux pl ane either (see, e.g., figure 4 in [Rodriguez et al.l 
l2011bl and references therein). 

We therefore made a rough estimate of the source 
flux in each observation by calculating the average-per- 
observation energy spectrum. In the upper panel of 
Figure [T] we plot the 2-50 keV unabsorbed flux. In 
IGR J17091-3624 we observe the a, /3, i/, p and p classes 
when the flux is between ~2 and '-^3.7 10~^ erg cm~^ s~^ 
(the diskbb temperature kT vary between 1 and 2 keV 
while the power law index between 2 and 3). As- 
suming that IGR J17091-3624 is emitting at Eddington 
rates, we derived the mass of the compact object as a 
function of distance (Figure [5|) for a 2-50 keV flux of 
4 • 10^^ erg cm~^ s~ . Given that the correction fac- 
tor linking the 2-50 keV and the bolometric flux is not 
exactly known, we plot curves for factors of 1, 2 and 3. 
Figure [S] implies that if IGR J 1709 1-3624 emits at Ed- 
dington, then either it harbors the lowest mass black hole 
known today (< 3Mq for distances lower than 17 kpc), 
or, it is very distant. Such a large distance, together with 
its b ~2.2° Galactic latitude, would imply a significant, 
but not necessarily implausible, altitude above the disk 
(e.g., b --2.8° and d > 25kpc for the BHC GS 1354-64, 
ICasares et al| [20091. 

Clearly, constraining the distance to IGR J17091-3624 
is fundamental for understanding the physical processes 
that govern the X-ray variability in IGR J17091-3624, 
and hence those in GRS 1915+105. The interstellar 
absorption in the direc tion of I GR J17091-362 4 is high 
(~ 1 X 10^^ cm~^, see iRodriguez et al]l2011bl ). imply- 
ing that optical studies when the b inary system is in 
quiescence (e.g., ICasares et al.ll200l ) might be challeng- 
ing. However, accurate parallax distance estimates like 
those reported for the bla ck hole X-ray binary V404 Cyg 
()Miller-Jones et al.|[2009l) could still be possible. 
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